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Исследуется связь между релятивистской термодинамикой ранней Вселенной с метрикой Логунова и гравитационным аналогом 
статистической механики. Обнаружено, что обратимое во времени уравнение Лиувилля может иметь частное решение с нару- 
шенной симметрией во времени. 


Существует мнение, что описание необратимых 
процессов на языке индивидуальных траекторий и 
описание на вероятностном языке не являются эк- 
вивалентными [1-3]. Для того чтобы вычислить 
энтропию, рассмотрим с точки зрения релятивист- 
ской теории гравитации Логунова связь между тер- 
модинамикой и гравитационным аналогом стати- 
стической механики. Одновременно это позволит 
исследовать проблему обратимости эволюции ран- 
ней Вселенной. 


Г равитационный аналог уравнения Лиувилля 


Сначала рассмотрим случай (ІІ 0 < 0), который 
при квантовом рассмотрении соответствует дис- 
кретному спектру. Гамильтониан гравитационного 
атома ( С 0 <0) имеет вид [4] 


Н(а,р а ) 
Ал а 2 

здесь V = — - — , а 
Логунова 


- -ѴрА, 
2М:а 1 01 


( 1 ) 


- масштабный фактор метрики 


сіст 2 = а 6 (х°)(сіх 0 ) 2 -а 2 (х°)х 


х[с!г 2 + г 2 ($т 2 (9)с1у/ 2 + <і& 2 )]. 


Ц)<0 - постоянная составляющая скалярного поля, 
р=—Мр(іа'— сопряженный координате а импульс, 

а' = ^-,сІ( = а і сІх 0 - собственное время. При этом 

переменные а(х°) и р а (х°) изменяются со временем 
х° согласно уравнениям Гамильтона 

да _ д г ГГ гг, д Р° - д 


дх° др а [ ^ Щ ’ дх° да 


= --~г~[\І8т Н], (2) 


здесь для метрики Логунова л/&» УР а " 

внений (2) следует, что эволюция фазовых точек во 
времени х° протекает в соответствии с теоремой 
Лиувилля. 

Для гамильтониана (1) ур. (2) можно предста- 
вить в виде 

да(і) _ дЙ др а ( 1 ) _ дЙ „) 

ді д р а ’ ді да 

здесь Н(а,р а ) = -^ Т - + Ѵр 0 \. 


Ур. (3) сводятся к одному дифференциальному 
уравнению второго порядка 

(— ) + 2(— ) 2 +8яС|С/ 0 | = 0, (4) 

а а 11 


полученному в работе [5]. Решение этого уравне- 
ния имеет вид 


а = а 0 сок 1/3 (ѵ I) 


и для такого решения 



Будем предполагать, что состояние может быть 
приписано гравитационному атому лишь с некото- 
рой вероятностью, т.е. по аналогии с квантовой 
теорией будем использовать вероятностную трак- 
товку динамических процессов одной частицы, а не 
ансамбля ее копий. Для этого без использования 
принципа сокращенного описания [3] введем 
функцию распределения р(а,р а ,(), которая пропор- 
циональна плотности вероятности расположения 
гравитационного атома в гравитационном фазовом 
пространстве ( а,р а ). 

Учитывая уравнения Гамильтона (2) и (3), не- 
трудно получить космологическое обобщение 
классического уравнения Лиувилля, описывающе- 
го эволюцию плотности р(а,р а ,1) в гравитационном 
фазовом пространстве (а,р а ) 


др у .д Н д р 

—=Ьр - 3 — , 

ді да др а 


( 5 ) 


у дН д дН д 

где ь = — — — — - оператор Лиувилля эф- 

да о р а о р а да 

фективного гамильтониана 0 . При этом, как и в 
классическом случае, функция распределения пос- 


тоянна вдоль фазовой траектории 


сі р 
сіі 


= 0 


Если функция р(а,р а ,і) описывает эволюцию, то в 
классической теории считается, что обращенная эво- 
люция описывается функцией р(а,-р„-1). Так как га- 
мильтониан Я обладает свойством Н(а,р а )=Н(а-р а ), 
то функция р(а,-р а ,-і) удовлетворяет ур. (5). Поэто- 
му это уравнение инвариантно относительно опера- 
ции обращения времени (а^а,р а ->-р а ,1->-і). Одна- 
ко, в соответствии со вторым началом термодинами- 
ки часто реальные системы не обладают симметрией 
по отношению к обращению времени, что противо- 
речит обратимости ур. (5) так же как, и других основ- 
ных уравнений классической и квантовой теории. 
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Тем не менее, существует связь между обратимой ми- 
кроскопической динамикой и необратимой макро- 
скопической эволюцией [1-3]. 

Для гамильтониана (1) решение ур. (5) имеет вид 


р(сі,р а ,і) = Р 


ѵі - агсІ§ 


н й му 


я. 


( 6 ) 


где Н„ = 


8я :ОЮ, 


°\н. = 4^н. 


Если предположить, что 


р(а, р а ,1) = р(1)' Р 
то тогда 

Р(а,р а ,1) = р(Н,)-ех р 
либо 


-агсТц 


Н 0 М]а 2 


,Я, 



с 

- 

у 

У 

ѵ і - агс1§ 

Ра 

Н„М 2 а 2 


_ 

V 

|_ ° р 2 

)_ 


,(7) 


р{а,р а ,() = рШ,) + 


+ ехр 


ѵ I - агсіа 


Н 0 М;а 2 


р(а,р а ,1) = 


= Р(Я.) 




( 

- 

\ 


1 + ехр 

У 

ѵі- агсІ§ 

V 

Ра 



Н 0 М 2 р а 2 _ 

/ 



( 8 ) 


где произвольная функция гамильтониана р(Щ со- 
ответствует равновесному состоянию, у = -8§п(г), 

-1, КО 

5§п(г) = <0, і = 0 ; величина у может быть и константой. 
1, 1>0 

Непосредственной подстановкой нетрудно убе- 
диться, что при ІІ 0 > 0 формальное решение гравита- 
ционного аналога уравнения Лиувилля (5) имеет вид 
р(а,р а ,1) = р(Н + ) + 


+ехр 



( 

- 

\ 

У 

ѵ 1 + Агііі 

Ра 


Н 0 М;а 2 _ 



_ 

V 

/_ 


( 9 ) 


либо 


Р(а,р а ,1) = 


= Р(Я + ) 



г 

г 

- 

\ 
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1+ехр 
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ѵ 1+ АгіЬ 

Ра 
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Н п М~ а ' 
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прошлым и будущим, и, что время является пара- 
метром удаленности от равновесного состояния. 
Поэтому статистическое описание на основе функ- 
ции распределения при у=-8§п(/) включает в себя 
приближение к равновесному состоянию и обрати- 
мость. Если же у=і, то решение обратимо, но не 
описывает приближение к равновесию. 

Когда же величина у является вещественной 
константой, то тогда р(а,р ю Г)фр(а,-р а ,-(). При 
этом одна из функций распределения стремится к 
равновесной функции р(Щ, а другая удаляется от 
нее, так что в этом случае прошлое и будущее раз- 
личимы. Таким образом, из инвариантности фун- 
даментальных законов относительно обращения 
времени не следует, что прошлое и будущее играют 
одинаковую роль. 

Будем считать реальным то распределение, ко- 
торое соответствует равновесному состоянию, до- 
стигаемому в будущем (т.е. при Т->со). Для этого в 
выражении (8) следует положить у=— 1. Описание 
на основе такой функции распределения включает 
в себя релаксацию к равновесному состоянию и 
необратимость. 


На фазовой же траектории, заданной в параме- 
трическом виде й=а 0 со8 1/3 (ѵО, р=-М 2 аа ' аргумент 
функции распределения 

Г \ 


I - агсГц 


у Н 0 М 2 р а 2 ; 


= 0, 


так что в этом случае функция распределения сов- 
падает с равновесной функцией 

Таким образом, обратимое во времени уравне- 
ние Лиувилля может иметь частное решение с нару- 
шенной симметрией во времени даже для индивиду- 
альной эффективной частицы Планка, обладающей 
массой М р ^ 10 5 г. Замечательно, что необратимость 
возникает в ситуации с одной степенью свободы. 
Это опровергает интерпретацию необратимости, со- 
гласно которой необратимость - не фундаменталь- 
ный закон природы, а следствие приближенного 
макроскопического характера наблюдений. 

Однако в обоих случаях при заданных началь- 
ных условиях можно с определенностью предска- 
зать будущее, или восстановить прошлое значение 
функции распределения, а значит и состояние си- 
стемы, т.е. детерминизм сохраняется. Если же 
предположить, что начальные данные являются 
случайными величинами, то можно приспособить 
классическую и квантовую теорию для описания 
таких процессов, для которых характерен выбор 
потенциальных возможностей. 


здесь Н + = 

Нетрудно видеть, что функция (8) и обращен- 
ная по отношению к ней функция р(а,-р а ,-1) рав- 
ны и в фиксированной точке фазового пространства 
за время 1>ѵ ' они стремятся к равновесной функ- 
ции р(Н.) Это означает эквивалентность между 


2 М;а 


-+ѴѴ п 


Энтропия ранней Вселенной 

Равновесную функцию распределения можно 
найти из принципа максимума информационной 
энтропии [3] 

5 ш (р) = -к ь \р(Н,)Щр(Н.У) ^ Ла сІр а , 
при дополнительных условиях 
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Я. = 


8 оо Р(Я,)ф я =1, 

—00 

00 

&Ю Л а | Я. р(Я.)ф я , 


вычисленных на фазовой плоскости. Численный 
расчет по области, ограниченной фазовой траекто- 
рией, совпадает с результатом по фазовой плоско- 
сти с точностью порядка единицы. Найдем экстре- 
мум функционала 

$(Р) = $ы~ а \р( н *)\[8^ Лси1р а - 
-Р\н,^<1ас1р а 

с множителями Лагранжа а и /3. Приравнивая ну- 
лю вариацию функционала, получим 

5§ = -^[1п(р) + рН,+а+1]8р^8^с1ас1р а =0. 

Т.к. 5р произвольна, то экстремальная функция 
распределения имеет вид канонического распреде- 
ления Гиббса 


Р(Я.) = ех р 


5(т) = ~|| Р\п (Р^сіа с}р а 

м 


ЗА. 

71 


х -1 
Хх 


{ х ъ йх | Р(х,у,х) Р{х,у,х)Ау + 8 0 

і_ Аѵ-і 

— Д 


где 


А(х,у,-г) = 


1п 


Ш + Ы2 ехр 


Г /НК + 


+агйі 


Л) 


+Ах~ 


V 


Р(х,у,т) = —ехр 


Г 

Л ,2 1\ 

Ах 3 

— -X 3 

V 

І_ х V 


#1 + N2 ■ ехр |Ѵ(т) т +агйі 


М=- 


д/я 2 х 6 - 


[ 71 й •] 

(Ю) 

31 ?хФх Т ехо 

яг 'Вф) 

2 со 

1 \ МЛ 1 слр 

V / 


1 / з 2 .ѵ 6 — 1 

з ГГ 





тгЧя. 
2 со 


и для гамильтониана (1) удовлетворяет нормиро- 
вочному условию |р(Я.)Ѵйо"Фф 0 = 1; множитель 

Лагранжа Д выражен через параметр ^ = Д, 
я/2со=Р, и его можно отождествить с температурой 
Т = к ь - постоянная Больцмана; Й=—Н„ Для 

энтропии из (10) следует выражение 
* 5 ' = -^,{ р(Н,)Ыр(Н,))^АаАр Л к 
где для гамильтониана (1) интеграл 

со оо 

(я) = | ^ | Я, р(Я.)ф в 

—00 —00 

легко вычисляется в аналитическом виде и равен 
(Й,^ = так что равновесное значение энтропии 

каждого отдельного гравитационного атома равно 

(11) 


А2 = 


ЗА 

л: 


| х ъ Ах | 


ехр 


Ял: 3 1 х 3 +агйі X 

х I х 


сіу 


Ш 


х = И?- ■ а, у- 

2\[к ’ 


/А 2 х 6 -1 

= -х 2 1Ь(агсЬ(А X 3 )) = - Ѵ — , $ 0 = 0, 060, 

Ах 

соответствующей распределению (9), представлена 
на рис. 1. Для определенности выбрано значение 
параметра А=0,01. Графики построены с помощью 
компьютерной программы «МаДюаё 2001». С по- 
мощью этой программы при каждом значении вре- 
мени г вычислялся двойной интеграл по перемен- 
ным а и р а по области фазового пространства, огра- 
ниченной фазовой траекторией. 


и для А 0 гравитационных атомов 

$ Ч =ЯА. (12) 

і 3 н 

Полагая к ь Т{ 0) = — = — — - и выбирая а а =2НР, 
л 2л 

ІІ й =ХМ р , константу 7Ѵ 0 в равенстве (12), являющую- 
ся аналогом числа Авогадро, умноженного на чи- 
сло молей, можно представить в виде /Ѵ п =1/А, 
А - безразмерная константа самодействия скаляр- 
ного поля. 

Нетрудно видеть, что при Ц > 0 функция распре- 
деления Гиббса (10) нормируема по области, огра- 
ниченной фазовой траекторией, а не по всей фазо- 
вой плоскости. Эволюция энтропии 



Рис. 1. Зависимость безразмерной энтропии на одну части 
цу от времени г -ѵі при у=~здп( I), и 0 >0 
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Отметим, что в этом случае для приближенного 
описания эволюции энтропии можно использовать 
удовлетворяющее принципу экстремума информа- 
ционной энтропии [3] сингулярное решение ура- 
внения (5) 

р = Ж- 1 5(Й + -Ё), 

где ]Ѵ=\\8(Й + -Ё)л[§^сІа(1р а , 5(х) - дельта-функция. 
Дело в том, что по аналогии с принципом локально- 
го равновесия в пространстве [6] естественно пред- 
положить, что в каждый момент времени существу- 
ет состояние равновесия. Проведенное выше рас- 
смотрение показывает, что в произвольный, но фик- 


сированный момент времени I объем V = 


4л а 


число частиц тѴ„ = 


4л аі 


и энергия Е являются 


0 Зк ь Т(0) 

заданными параметрами. Поэтому в каждый мо- 
мент времени энтропию можно определить по фор- 
муле Больцмана 

8(Ё,Ы,Ѵ)=]п(ЩЁ,Ы,Ѵ)), 

где 1Ѵ=\\8(Й + -Ё)^сІа(Ір а - статистический вес, 
задающий число динамических состояний внутри 
слоя, определяемого неравенствами Ё<Й+<Ё+АЕ', 




РІ 


-Ѵр» 


2 М\а 

решения (1) нетрудно получить 

ЩЁ,И,Ѵ) = — АгкЬ 
2(о 


>\/Ім - тг [5]. Для 


ѵ 1 

[Щощ, 

ы 

Ё 


здесь Ѵ 0 = 


4л а: 


. Подставляя в последнее выраже- 


4л а„ 


ние V = ° созЬ(уг), получим выражение, при- 

ближенно (г<1) описывающее эволюцию энтро- 
пии Больцмана ранней Вселенной, 


5(т) = 1п 


Аг 8 І 1 


созЬ(г) 


кТ(0)М 0 




где (т=ѵІ), а величина л -1 интерпретируется как воз- 
раст Вселенной. При этом для времен г < 1 ^ » 0. 


Эволюция безразмерной энтропии Больцмана 
представлена на рис. 2. 



Рис. 2. Зависимость безразмерной энтропии Больцмана от 
времени т 


Из рис. 1 и 2 видно, что эволюция системы име- 
ет смысл только при ?>0 и энтропия из-за неустой- 
чивости решения (1) может возрастать беспредель- 
но с малым темпом возрастания. Отсутствие мак- 
симума энтропии обусловлено возрастанием объе- 
ма V. Ограниченность энтропии может обеспечить 
скалярный потенциал вида 


Щ<Р) 


Г^о. 

1 °, *>Ё 


Момент времени і т можно определить из усло- 
вия 5(4)=1. В этом случае 4»1,8ѵ 1 , а параметр са- 
модействия скалярного поля А=1 и, следовательно, 
й=1. Из рис. 1 и 2 видно, что время релаксации не- 
стационарного состояния /^2ѵЛ 

Эволюция безразмерной энтропии отдельной 
частицы для распределения (8) представлена на 
третьем и четвертом рисунках, где т=ѵ/, Я= 1,18, 
а 0 =2НЁ- Значение константы теории А выбрано из 
условия 5(0)=0, являющегося аналогом третьего 
закона термодинамики. Выбор начала отсчета вре- 
мени ? 0 =0 обусловлен экстремумом решения ур. (4). 



Рис. 3. Зависимость безразмерной энтропии на одну части- 
цу от времени х=ѵі при у=~1 


Из рис. 3 видно, что энтропия возрастает как при 
т<0, так и при т>0, но в первом случае 5”<0, а во вто- 
ром случае 5>0. Поэтому естественно предположить, 
что эволюция системы имеет смысл только при />0. 



X 

Рис. 4. Зависимость безразмерной энтропии на одну части- 
цу от времени т=ѵі при у=~здп(1), ІІ 0 <О 


В заключение отметим, что для потенциала 

-Ѵ 0 ,0 <(<(, 

о, і > ( 2 
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рассмотренные в статье случаи ((/ 0 > 0) и (1/ 0 < 0) 
можно рассматривать как этапы единого процесса 
эволюции ранней Вселенной. 

Интересно отметить, что на линии постоянного 
уровня Н,(а,р)=Е распределение (10) совпадает с 
вероятностью рождения Вселенной, найденной в 
работе [7]. 

Из (10) видно, что каноническое распределение 
Гиббса не является аналитичным по константе взаи- 
модействия со. Поэтому оно не может быть получе- 
но на основе теории возмущений. Функция распре- 
деления (10) больше там, где |7/ 0 | имеет большие зна- 
чения. Если предположить, что | Щ является случай- 
ной величиной, то тогда функция распределения 
станет случайной функцией. По определению слу- 
чайной функцией называют функцию неслучайного 
аргумента і, которая при каждом фиксированном 
значении аргумента является случайной величиной. 
Такие функции позволят описывать индетермини- 
стические процессы, характерные для живой приро- 
ды. Такой подход осуществлен в работе [8]. 

Отметим, что в соответствии с ур. (5) всегда пол- 
ная производная ^ = 0, а частная производная мо- 
жет отличаться от нуля. Поэтому для любого распре- 
деления, удовлетворяющего гравитационному ана- 
логу ур. Лиувилля (5), полная производная инфор- 
мационной энтропии ^ = 0, а * 0. Частная 
производная определяет эволюцию энтропии в 
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фиксированной точке фазового пространства, что 

может иметь место в процессе рождения Вселенной. 

Заключение 

1. Получено космологическое обобщение класси- 
ческого уравнения Лиувилля, и найдено точное 
его решение, описывающее релаксацию состоя- 
ния гравитационного атома к его равновесному 
состоянию. Обнаружено, что детерминистиче- 
ское уравнение может иметь частное решение с 
нарушенной симметрией во времени. Найдены 
решения, описывающие как обратимые, так и 
необратимые процессы. Все это может служить 
динамическим обоснованием необратимости. 

2. Получено уравнение состояния гравитацион- 
ного атома. Найден гравитационный аналог 
первого закона термодинамики. 

3 . Вычислена в аналитическом виде энтропия рав- 
новесного и численно нестационарного состоя- 
ния гравитационного атома, безразмерная эн- 
тропия которого возрастает до тех пор, пока в 
состоянии равновесия не достигнет максималь- 
ного значения, равного единице. 

4. Рождение обычного вещества во Вселенной, ли- 
бо мини-вселенных - событие, связанное с на- 
рушающей симметрию во времени неустойчиво- 
стью гравитационного атома относительно 
спонтанных переходов. Гравитационные атомы 
могут быть той средой, которая породила обыч- 
ное вещество Вселенной, либо мини-вселенные. 
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